■  ^  ■  Defence  Research  and  Recherche  et  ddveloppement 
m  ~  u  Development  Canada  pour  la  defense  Canada 


Micro-Doppler  Analysis  of  Rotating 
Target  in  SAR 


T.  Thayaparan,  S.  Abrol  and  S.  Qian 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


Defence  R&D  Canada  -  Ottawa 

TECHNICAL  MEMORANDUM 
DRDC  Ottawa  TM  1005-204 
December  2005 


Canada 


Micro-Doppler  Analysis  of  Rotating  Target  in 
SAR 


T.  Thayaparan 

Defence  R&D  Canada  -  Ottawa 
S.  Abrol 

University  of  Ottawa 
S.  Qian 

National  Instruments 


Defence  R&D  Canada  -  Ottawa 

Technical  Memorandum 
DRDC  Ottawa  TM  2005-204 
December  2005 


©  Her  Majesty  the  Queen  as  represented  by  the  Minister  of  National  Defence,  2005 
©  Sa  majeste  la  reine,  representee  par  le  ministre  de  la  Defense  nationale,  2005 


Abstract 


Rotating  targets  cause  phase  modulation  of  the  azimuthal  phase  history  of  a  SAR 
system.  The  phase  modulation  may  be  seen  as  a  time-dependent  micro-Doppler  (m- 
D)  frequency.  Due  to  their  superior  resolution  potential,  it  is  useful  to  analyze  such 
signals  with  time-frequency  analysis  methods.  This  report  presents  two  approches 
for  extracting  m-D  features  from  SAR  images.  In  order  to  extract  m-D  features 
from  SAR  images,  the  time  domain  radar  return  is  decomposed  in  two  separate 
ways.  One  is  based  on  wavelet  decomposition  in  which  the  returned  signal  is  de¬ 
composed  into  a  set  of  components  that  are  represented  at  different  wavelet  scales. 
The  components  are  then  reconstructed  by  applying  the  inverse  wavelet  transform. 
This  wavelet  approach  has  been  used  in  our  previous  m-D  analysis  work  for  a  ISAR 
system,  and  it  is  presented  here  in  the  extraction  of  m-D  features  for  a  SAR  sys¬ 
tem.  The  second  approach  is  based  on  adaptive  chirplet  decompostion.  This  new 
approach  is  introduced  as  an  alternative  to  the  wavelet  approach  of  decomposing 
the  SAR  radar  return.  The  results  from  the  wavelet  and  adaptive  chirplet  decom¬ 
postion  procedures  are  compared,  and  the  chirplet-based  approach  establishes  itself 
as  a  viable  alternative.  The  chirplet  based  method  of  m-D  extraction  has  been  suc¬ 
cessfully  applied  to  SAR  data  scene  collected  by  the  US  Navy  APY-6  radar.  The 
rotating  frequency  of  the  antenna  was  calculated  using  this  method  and  the  results 
agree  well  with  available  ground  truth  data. 
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Resume 


Les  cibles  rotatives  causent  la  modulation  de  la  phase  azimutale  d’un  systeme  SAR. 
La  modulation  de  phase  peut  etre  consideree  comme  une  frequence  micro-Doppler 
(m-D)  variable  dans  le  temps.  Vu  leur  potentiel  de  resolution  superieure,  il  est 
utile  d’analyser  ces  signaux  au  moyen  de  methodes  d’analyse  temps-frequence.  Le 
present  rapport  decrit  deux  methodes  d’extraction  des  caracteristiques  m-D  des  im¬ 
ages  SAR.  En  vue  d’extraire  les  caracteristiques  m-D  des  images  SAR,  l’echo  radar 
dans  le  domaine  temps  est  decompose  de  deux  fayons  distinctes.  L’une  est  basee  sur 
la  decomposition  de  l’ondelette  dans  laquelle  l’echo  est  decompose  en  un  ensemble 
de  composantes  qui  sont  representees  sur  differentes  echelles  d’ondelette.  Les  com- 
posantes  sont  ensuite  reconstitutes  par  l’application  de  la  transformee  d’ondelette 
inverse.  Cette  methode  a  ondelette  a  ete  utilisee  dans  notre  travail  d’analyse  m-D 
anterieur  pour  un  systeme  ISAR,  et  elle  est  presentee  ici  dans  I’extraction  des  ca¬ 
racteristiques  d’un  systeme  SAR.  La  seconde  methode  est  basee  sur  la  decompos¬ 
ition  de  Chirplet  adaptatif.  Cette  nouvelle  methode  est  introduite  comme  rechange 
de  la  methode  a  ondelette  pour  decomposer  1’echo  du  radar  SAR.  Les  resultats  des 
procedures  de  decomposition  d’ondelette  et  de  Chirplet  adaptatif  sont  compares,  et 
la  methode  basee  sur  Chirplet  s’avere  une  methode  de  rechange  viable.  La  methode 
basee  sur  Chirplet  pour  extraire  les  caracteristiques  m  D  a  ete  appliquee  avec  suc- 
ces  a  la  collecte  de  donnees  de  scene  SAR  au  moyen  du  radar  APY-6  de  la  mar¬ 
ine  americaine.  La  frequence  de  rotation  de  l’antenne  a  ete  calculee  au  moyen  de 
cette  methode,  et  les  resultats  concordent  bien  avec  les  donnees  de  realite  de  terrain 
disponibles. 
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Executive  summary 


SAR  images  are  a  high-resolution  map  of  surface  target  areas  and  terrain  in  the 
range  and  cross-range  dimension.  If  there  are  moving  targets  in  the  scene,  SAR 
cannot  simultaneously  produce  clear  images  of  both  these  stationary  targets  and 
moving  targets.  Usually,  moving  targets  appear  as  defocussed  and  spatially  dis¬ 
placed  objects  superimposed  on  the  SAR  map.  Moving  targets  affect  the  phase 
history  of  a  SAR  collection  in  a  manner  dependent  on  the  particular  motion  of  the 
target.  Some  targets  contain  parts  that  move  relative  to  the  target  itself  Examples 
are  rotating/vibrating  parts  such  as  wheels  or  engines.  A  special  case  is  rotating 
point  targets  such  as  a  rotating  antenna.  The  resulting  class  of  motions  may  aid 
target  characterization  and  recognition.  Rotations  and  vibrations  can  be  observed 
by  radar  when  the  conditions  are  right.  The  phenomenon,  as  observed  by  radar,  is 
termed  micro-Doppler  (m-D).  In  addition  to  the  m-D,  there  may  be  a  Doppler  shift 
corresponding  to  the  target  body  motion. 

In  this  report,  we  examine  the  case  of  a  rotating  antenna  on  the  ground.  Such  a  ro¬ 
tating  target  causes  sinusoidal  phase  modulation  of  the  azimuth  phase  history  of  the 
SAR  collection.  The  phase  modulation  may  equivalently  be  seen  as  a  time-varying 
Doppler  frequency.  Time-frequency  methods  are  useful  tools  in  the  analysis  of  such 
signals  as  they  use  time  integration  while  still  allowing  for  non-stationary  signals. 
Also,  multicomponent  signals  may  be  analyzed.  Therefore,  a  method  for  extracting 
target  motion  parameters  is  needed  so  that  motion  parameters  may  be  estimated  and 
targets  can  be  identified. 

This  report  presents  two  approches  for  extracting  m-D  features  from  SAR  images. 
In  order  to  extract  m-D  features  from  SAR  images,  the  time  domain  radar  return  is 
decomposed  in  two  separate  ways.  One  is  based  on  wavelet  decomposition  in  which 
the  returned  signal  is  decomposed  into  a  set  of  components  that  are  represented  at 
different  wavelet  scales.  The  components  are  then  reconstructed  by  applying  the 
inverse  wavelet  transform.  This  wavelet  approach  has  been  used  in  our  previous  m- 
D  analysis  work  for  a  ISAR  system,  and  it  is  presented  here  in  the  extraction  of  m-D 
features  for  a  SAR  system.  The  second  approach  is  based  on  the  adaptive  chirplet 
decompostion.  This  new  approach  is  introduced  as  an  alternative  to  the  wavelet 
approach  of  decomposing  the  SAR  radar  return.  The  results  from  the  wavelet  and 
adaptive  chirplet  decompostion  procedures  are  compared,  and  the  chirplet  based 
approach  establishes  itself  to  be  a  viable  alternative.  The  chirplet-based  method 
of  m-D  extraction  has  been  successfully  applied  to  SAR  data  scene  collected  by 
the  US  Navy  APY-6  radar.  The  rotating  frequency  of  the  antenna  was  calculated 
using  this  method  and  the  results  agree  well  with  available  ground  truth.  Such 
parameters  may  be  useful  for  describing  the  target,  or  serve  as  inputs  to  automatic 
target  recognition  algorithms. 
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Sommaire 


Les  images  SAR  sont  une  carte  a  haute  resolution  des  aires  de  cibles  de  surface  et 
du  terrain  sur  le  plan  de  la  distance  et  de  la  largeur.  Si  la  scene  comprend  des  cibles 
mobiles,  le  SAR  ne  peut  pas  produire  simultanement  des  images  claires  a  la  fois  des 
cibles  fixes  et  des  cibles  mobiles.  D’habitude,  les  cibles  mobiles  se  presentent  sous 
la  forme  d’objets  defocalises  et  deplaces  dans  l’espace,  superposes  a  la  carte  SAR. 
Les  cibles  mobiles  influent  sur  la  discordance  de  phase  d’une  collecte  SAR  d’une 
fagon  dependant  du  mouvement  particulier  de  la  cible.  Certaines  cibles  contiennent 
des  parties  qui  se  deplacent  par  rapport  a  la  cible  elle-meme.  Des  exemples  sont  des 
parties  rotatives  ou  vibratoires  telles  que  les  roues  ou  les  moteurs.  Un  cas  special 
est  constitue  par  les  cibles  a  point  de  rotation,  comme  les  antennes  rotatives.  La 
classe  de  mouvements  qui  en  resulte  peut  aider  la  caracterisation  et  la  reconnais¬ 
sance  des  cibles.  Le  radar  permet  d’observer  les  rotations  et  les  vibrations  lorsque 
les  conditions  sont  appropriees.  Ce  phenomene,  observe  au  radar,  est  appele  effet 
micro-Doppler  (m-D).  En  plus  de  l’effet  m-D,  il  peut  y  avoir  un  decalage  Doppler 
correspondant  au  mouvement  de  l’objet  cible. 

Dans  le  present  rapport,  nous  examinons  le  cas  d’une  antenne  rotative  au  sol.  Une 
telle  cible  rotative  cause  la  modulation  sinusoidale  de  la  phase  azimutale  de  la  col¬ 
lecte  SAR.  La  modulation  de  phase  peut  egalement  etre  consideree  comme  une 
frequence  Doppler  variable  dans  le  temps.  Les  methodes  temps-frequence  sont  des 
outils  utiles  pour  analyser  ces  signaux,  car  elles  font  appel  a  1’integration  temporelle 
tout  en  admettant  des  signaux  d’objets  mobiles.  De  plus,  elles  permettent  d’analyser 
des  signaux  a  composantes  multiples.  Par  consequent,  il  faut  une  methode  pour  ex- 
traire  les  parametres  de  mouvement  des  cibles  pourqu’on  puisse  estimer  les  valeurs 
des  parametres  de  mouvement  et  identifier  les  cibles. 

Le  present  rapport  decrit  deux  methodes  d’extraction  des  caracteristiques  m-D  des 
images  SAR.  En  vue  d’extraire  les  caracteristiques  m-D  des  images  SAR,  I’echo 
radar  dans  le  domaine  temps  est  decompose  de  deux  fagons  distinctes.  L’une  est 
basee  sur  la  decomposition  de  I’ondelette  dans  laquelle  I’echo  est  decompose  en 
un  ensemble  de  composantes  qui  sont  representees  sur  differentes  echelles  d’onde- 
lette.  Les  composantes  sont  ensuite  reconstitutes  par  1’application  de  la  transformee 
d’ondelette  inverse.  Cette  methode  a  ondelette  a  ete  utilisee  dans  notre  travail  d’ana- 
lyse  m-D  anterieur  pour  un  systeme  ISAR,  et  elle  est  presentee  ici  dans  I’extraction 
des  caracteristiques  d’un  systeme  SAR.  La  seconde  methode  est  basee  sur  la  de¬ 
composition  d’un  Chirplet  adaptatif.  Cette  nouvelle  methode  est  introduite  comme 
rechange  de  la  methode  a  ondelette  pour  decomposer  l’echo  du  radar  SAR.  Les 
resultats  des  procedures  de  decomposition  d’ondelette  et  de  Chirplet  adaptatif  sont 
compares,  et  la  methode  basee  sur  Chirplet  s’avere  une  methode  de  rechange  viable. 
La  methode  basee  sur  Chirplet  pour  extraire  les  caracteristiques  m  D  a  ete  appli- 
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quee  avec  succes  a  la  collecte  de  donnees  de  scene  SAR  au  moyen  du  radar  APY-6 
de  la  marine  americaine.  La  frequence  de  rotation  de  l’antenne  a  ete  calculee  au 
moyen  de  cette  methode,  et  les  resultats  concordent  bien  avec  les  donnees  de  realite 
de  terrain  disponibles.  Ces  parametres  peuvent  etre  utiles  pour  decrire  la  cible  ou 
servir  de  valeurs  d’entree  des  algorithmes  de  reconnaissance  automatique  de  cibles. 
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1  Introduction 


1.1  Synthetic  Aperture  Radar 

The  concept  of  radar  is  actually  quite  simple.  Microwave  radiation  is  transmit¬ 
ted  to  a  target  location,  and  the  reflected  energy  is  received  and  measured  by  an 
antenna/receiver  system.  This  idea,  however  simple,  is  applied  in  a  tremendous 
variety  of  rather  complex  applications.  One  field  of  application  is  broad-area  ima¬ 
ging.  In  such  a  system,  the  resolution  is  affected  by  the  size  of  the  aperture,  (i.e., 
the  area  used  to  receive  signals);  a  larger  aperture  results  in  a  higher  resolution.  In 
other  words,  radar  performance  is  inhibited  by  the  size  of  the  emitting  antenna.  A 
much  larger  antenna  than  can  be  practically  employed  is  needed  for  an  acceptable 
resolution.  If  the  radar  is  equipped  to  a  moving  platform,  such  as  a  satellite  or  an 
aircraft,  then  it  is  possible  to  combine  reflected  signals  from  along  the  flight  path 
to  simulate  (synthesize)  a  very  long  antenna.  The  radar  transmits  a  wide  beam  of 
electromagnetic  signals  at  a  series  of  points  along  the  flight  path.  Since  the  signals 
are  coherent  and  always  in  phase,  the  radar  return  is  collected  and  then  combined 
in  order  to  create  high  resolution  images  just  as  they  would  have  been  made  if  a 
very  large  antenna  had  been  employed.  In  this  way,  a  long  aperture  is  artificially 
created.  This  is  known  as  Synthetic  Aperture  Radar  (SAR).  Here,  it  is  important  to 
note  that  the  radar  is  moving  relative  to  the  target.  This  is  quite  contrary  to  ISAR 
(Inverse  Synthetic  Aperture  Radar)  where  an  image  is  generated  from  reflected  sig¬ 
nals  collected  by  a  stationary  radar  from  a  moving  target  [1-7]. 

SAR  systems  usually  produce  two-dimensional  images.  One  dimension  is  called 
the  range  and  is  obtained  by  precisely  measuring  the  time  from  transmission  of  a 
signal  to  its  return.  The  second  dimension  is  known  as  the  azimuth  and  is  perpen¬ 
dicular  to  the  range  dimension.  The  azimuthal  resolution  is  obtained  by  processing 
the  Doppler  phase  of  the  radar  return.  What  differentiates  SAR  from  other  radars  is 
its  ability  to  produce  fine  azimuth  resolution  [1-7]. 

SAR  is  a  well  established  and  useful  technique  of  acquiring  high-resolution  images 
of  an  area  of  interest  from  airborne  or  space  sensors  [1-2].  A  particular  class  of 
targets  that  pose  a  difficult  challenge  for  target  recognition  is  moving  targets.  If 
there  are  moving  targets  in  the  scene,  SAR  cannot  simultaneously  produce  clear 
images  of  both  of  these  stationary  targets  and  moving  targets  [3-7].  Usually,  moving 
targets  appear  as  defocused  and  spatially  displaced  objects  superimposed  on  the 
SAR  map.  Moving  targets  cause  phase  modulation  of  the  azimuth  phase  history  of 
a  SAR  collection.  The  phase  modulation  can  be  seen  as  a  time-varying  Doppler 
frequency.  It  is  obvious  that  in  some  scenarios  it  is  the  moving  objects  that  are  of 
interest.  Therefore,  a  method  for  extracting  target  motion  parameters  is  needed  so 
that  motion  parameters  may  be  estimated  and  targets  can  be  identified. 
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1.2  Micro-Doppler  Phenomenon 


When  a  radar  transmits  an  electromagnetic  signal  to  a  target,  the  signal  interacts 
with  the  target  and  then  returns  to  the  radar.  Clearly,  any  changes  in  the  properties  of 
the  returned  signal  are  generally  due  to  its  interaction  with  the  target.  This  being  the 
case,  the  characteristics  of  the  target  in  question  are  reflected  in  the  returned  signal, 
and  it  follows  that  various  properties  of  the  relevant  target  can  be  extracted  from  the 
returned  signal.  For  example,  when  the  target  is  moving,  the  carrier  frequency  of 
the  returned  signal  will  be  shifted  due  to  the  Doppler  effect.  This  Doppler  frequency 
shift  can  then  be  used  to  determine  the  radial  velocity  of  the  moving  target. 

It  frequently  occurs  that  a  target  or  some  structure  on  the  target  is  vibrating  or  rotat¬ 
ing  in  addition  to  the  target’s  translation.  These  vibrations  and  rotations  are  referred 
to  as  micro-motion  dynamics.  The  micro-motion  dynamics  of  a  target  generate  fre¬ 
quency  modulations  on  the  returned  signal  in  addition  to  the  Doppler  frequency 
shift  caused  by  the  target’s  translational  motion.  In  fact,  the  frequency  modula¬ 
tions  on  the  returned  signal  induced  by  the  micro-motions  produce  sidebands  about 
the  target’s  Doppler  frequency.  The  frequency  modulation  due  to  micro-motion 
dynamics  is  called  the  micro-Doppler  (m-D)  phenomenon  [7-13], 

While  the  Doppler  frequency  shift  created  by  the  translational  motion  of  the  target 
is  time-invariant  at  a  constant  velocity,  the  Doppler  frequency  shift  generated  by 
micro-motion  dynamics  is  a  time-varying  function  that  imposes  a  periodic  time- 
varying  modulation  onto  the  earner  frequency.  The  modulation  contains  harmonic 
frequencies  that  depend  on  the  carrier  frequency,  the  vibration  or  rotation  rate,  and 
the  angle  between  the  direction  of  vibration  and  the  direction  of  the  radar’s  incident 
wave. 

The  m-D  phenomenon  is  of  interest  because  information  regarding  a  target’s  micro¬ 
motion  dynamics  is  preserved  in  the  returned  radar  signal.  The  m-D  phenomenon 
is  commonly  encountered  in  radar  returned  signals  because  real  world  targets  (e.g., 
helicopter,  vehicle,  rotating  antenna)  are  usually  engaged  in  complicated  maneuvers 
that  incorporate  translation,  vibrational  and  rotational  motions.  The  mechanical 
vibration/rotation  of  a  target  or  part  of  the  target’s  structure  appears  in  a  radar  image 
as  smeared  features  across  the  cross-range  dimension  of  the  image  and  induces 
frequency  modulation  on  radar  returned  signals  [7-17],  On  the  other  hand,  m-D 
features  could  provide  distinctive  information  for  recognition  of  targets  of  interest. 
They  have  been  used  to  identify  the  natural  resonant  frequency  of  a  tractor- trailer 
truck  [14].  The  m-D  features  of  a  Jet  Engine  Modulation  (JEM)  lines  in  a  Mi- 
24  Hind-D  helicopter  are  also  used  to  estimate  the  turbine  rotation  rate  and  the 
number  of  turbine  blades  [15].  Also,  research  has  been  conducted  on  radar  Doppler 
signatures  in  the  area  of  human  gait  analysis  [7,9,1 1,13,  18-19], 
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It  is  reasonable  to  expect  that  the  m-D  features  representing  the  micro-motions  of  a 
target  can  be  extracted  from  the  returned  signal,  much  in  the  same  way  as  properties 
are  extracted  from  radar  returns  of  targets  undergoing  only  translational  motion. 
Since  different  targets  produce  different  micro-motions,  every  target  would  have  its 
own  ”m-D  signature”,  making  it  possible  to  distinguish  and  identify  targets  under 
consideration  based  on  the  additional  information  provided  by  the  m-D  features. 
Hence,  an  effective  method  is  needed  for  extracting  m-D  features  in  order  to  fully 
exploit  the  additional  and  unique  information  they  provide  [11], 

In  this  work,  we  extract  the  m-D  features  relating  to  a  rotating  antenna  in  a  SAR 
target  scenarios.  The  m-D  for  such  rotating  target  may  be  seen  as  a  sinusoidal 
phase  modulation  of  the  SAR  azimuth  phase  history.  Along  with  the  wavelet  de¬ 
composition  method  used  in  previous  works  of  m-D  analysis  [1 1-13],  we  introduce 
an  adaptive  chirplet  decomposition  method  to  extract  m-D  features.  The  goal  of  this 
report  is  to  find  a  method  of  extracting  m-D  features  in  SAR  and  detect  useful  tar¬ 
gets.  These  m-D  features  may  provide  additional  information  for  target  recognition 
complementary  to  existing  recognition  methods. 
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2  Method  of  Micro-Doppler  Feature 
Extraction 


Joint  time-frequency  analysis  is  the  base  of  most  of  the  existing  methods  used  to 
extract  m-D  features  [7-13],  The  procedure  of  wavelet  analysis  incorporated  with 
time-frequency  analysis  (specifically,  the  short-time  Fourier  transform)  has  been 
applied  in  our  previous  works  for  a  ISAR  system  [11-13].  Here,  the  approach 
is  used  in  a  SAR  target  scenarios.  Also,  we  introduce  a  new  procedure  for  m-D 
analysis,  namely,  the  adaptive  chirplet  decomposition  incorporated  with  the  time- 
frequency  distribution  series. 


2.1  Wavelet  Analysis 

A  very  useful  tool  in  many  recent  signal  processing  applications  such  as  detecting 
the  discontinuity  of  a  signal,  de-noising  a  signal,  and  image  compression,  the  wave¬ 
let  transform  results  in  an  efficient  representation  of  highly  non-stationary  signals 
[20-24].  The  main  motivation  for  applying  wavelet  analysis  in  the  extraction  of  m- 
D  features  is  that  the  micro-motion  dynamics  of  a  target  that  induce  m-D  features 
change  at  a  rate  much  faster  than  the  target  body  itself.  Wavelet  analysis  has  the 
capability  of  detecting  rapid  changes  of  a  signal  [20-24].  Therefore,  wavelet  ana¬ 
lysis  seems  ideal  for  the  extraction  of  m-D  features  and  the  wavelet  transform  can 
be  considered  a  powerful  tool  for  the  task. 

The  wavelet  analysis  is  actually  quite  complex  and  contains  a  large  number  of  steps. 
The  wavelet  transform  is  applied  by  using  a  tree  of  digital  filter  banks  based  upon 
the  multi-resolution  analysis  theory  [20-24].  In  order  to  effectively  extract  m-D 
features  from  the  returned  signal,  a  four  level  decomposition  tree  is  utilized  to  rep¬ 
resent  the  returned  radar  signal.  Using  the  wavelet  transform,  the  signal  is  first 
broken  down  into  its  constituent  parts  consisting  of  the  four  detail  levels  containing 
the  m-D  features  and  the  left  over  stationary  approximation  representing  the  body 
vibration.  The  inverse  wavelet  transform  is  then  used  to  reconstruct  the  decom¬ 
posed  wavelet  components.  During  the  reconstruction  process  of  the  detail  levels, 
only  the  wavelet  coefficients  that  are  related  to  the  m-D  features  of  the  signal  are 
used  while  the  other  coefficients  are  set  to  zero.  The  same  process  is  also  used 
in  the  reconstruction  of  the  stationary  part  of  the  target’s  body.  The  result  of  the 
wavelet  analysis  is  the  separation  of  the  target’s  micro-motion  dynamics  and  the 
extraction  of  the  m-D  features.  This  method  has  been  tested  previously  and  shown 
to  be  successful  with  other  data  sets  [11-15]. 

The  tree  of  digital  filter  banks  for  computing  the  discrete  wavelet  transform  is  given 
in  Figure  1  (this  is  a  four  level  decomposition  tree).  L  and  H  represent  pairs  of 
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Figure  1:  The  tree  of  filter  banks  for  computing  the  discrete  wavelet  transform. 


discrete  low-pass  and  high-pass  filters.  As  demonstrated  in  this  figure,  the  original 
signal  S  is  decomposed  into  its  constituent  parts  consisting  of  D\  ,  D2,  D3,  D4,  and 
A4.  In  other  words,  after  decomposition  S  =  A4  +  D\  +  D2.  +  D3  +  D4,  Now, 
the  wavelet  transform  has  decomposed  the  signal  into  five  parts;  the  four  detail 
levels  and  the  stationary  approximation.  Generally,  these  decomposed  parts,  by 
themselves,  do  not  represent  motion  dynamics  of  the  target  in  question.  Depending 
on  the  physical  quantity  that  one  wishes  to  estimate,  the  decomposed  parts  of  the 
signal  may  need  to  be  combined  in  various  ways  in  order  to  represent  the  m-D 
features  for  which  motion  parameters  are  required. 


2.2  Adaptive  Chirplet  Analysis 

Here,  an  alternative  to  wavelet  analysis,  the  so-called  adaptive  chirplet  analysis, 
is  proposed  in  the  extraction  of  m-D  features  from  a  radar  return.  A  chirplet  is  a 
windowed  portion  of  a  chirp.  A  chirp  being  a  function  very  similar  to  a  harmonic 
function  but  instead  of  having  a  fixed  period,  the  period  changes  with  time  along 
the  function.  The  chirplet  is  to  a  wavelet  as  the  chirp  is  to  a  wave,  as  shown  in 
Figure  2.  Many  events,  especially  those  found  in  nature,  can  be  modelled  as  the  su¬ 
perposition  of  short-lived  chirp  functions.  The  chirplet  transform  has  been  recently 
used  in  areas  such  as  multi-component  signal  detection,  instantaneous  frequency 
estimation,  time-varying  filter  design,  and  interference  excision  [9,  20,  25]. 


In  this  work,  the  adaptive  Gaussian  chirplet  decomposition  is  employed  to  extract 
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WAVE  WAVELET 


CHIRP  CHIRPLET 

Figure  2:  Chiplet  versus  Wavelet 


m-D  features.  The  Gaussian  chirplet  is  defined  as  [20,25] 

hk{t)  =  \/^exp  [  -3*  (t  -  tkf  +  j  (uk  {t  -  tk)  +  f  (t  -  tkf)  } 

ak  >  0  (0 

ik>Wk>0k  ^  ^ 

where  (tk)  uk)  indicates  the  time  and  frequency  center  of  the  linear  chirp  function, 
the  variance  ak  controls  the  width  of  the  chirp  function,  and  the  parameter  3 k  de¬ 
termines  the  rate  of  change  of  frequency.  When  the  chirp  rate  3k  is  equal  to  zero 
and  the  parameter  ak  approaches  zero,  the  linear  chirplet  function  hk(t)  reduces  to 
a  sinusoidal  signal. 

The  adaptive  spectrogram  for  the  Gaussian  chirplet  is  given  as  [20] 

2  -  —  (u  -  Ujk  -  3kt.)2  \  (2) 

otk  J 

The  algorithm  used  here  estimates  the  optimal  elementary  functions  in  the  adaptive 
signal  decomposition.  Converting  the  optimization  process  to  a  traditional  curve¬ 
fitting  problem  is  the  basic  idea  behind  this  algorithm.  The  relationship  between 


AS(t,uj)  =  2j2\Bk 


-ctk  ( t  -  tk) 
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the  estimated  parameters  and  testing  variables  is  given  by  the  function 


P  (&-ki  tki  k >  i  ^ k,rn  i  $k,in) 


(3) 


as  follows 


P  !^k'i  &k,mt  tk,mi  ^k,mj  Pk ,m) 

l-^-l 


(  - -  l  [(afc+«fc,m)2+(-9fc-3fc,m)2] 

x  =xp  I  (**  - 


X  exp 


A  m  *'"* - v?  (^fc  “  £&vm)  (^A*  — 


.  Qi--h Ocl-  riL 

x  exp  <;  -r- — : - —3-  5 


2  [(ttfr +£**... m)  +(^fc-^A-.r«)  ] 


TT  kifejp  ) 


(4) 


This  can  be  written  as 


P  t} c,  cj/c,  3k ;  q^q.  0) 

(^.,„)e-6tC7(u-‘)-t'fc-"‘]2 


where 


a  (ffc,»,)  = 


:  exp 


[(rtk+Qn.o)  +  (3*  —  ^-,o)  ] 


\J  V  4a*i0a*  y 

^  l+k un)  —  ^k  c*k+a k.o  ^/c,7ij 


(k-4,,,)2  \ 
2(nf  1+,Vo)  J 


(5) 


and 
6  = 


_ aA  +0^,0 _ 

2^(afc-fatio)  +  (#*  — #fc,o)  ] 


Since  both  sides  of  the  previous  equation  are  greater  than  zero,  we  have 


(6) 
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The  following  three  test  points  are  inserted  into  equation  5. 

Oik, 0,  ikfi,^k,-U0k  o) 

ttfc.o,  tkto,u>k,<hPkfi)  (7) 

This  results  in  the  linear  system  given  below 


(wfci0  -Wfc,  1 )  1  ) 

_  b 

I  (*k>hnk'0Jk'0,. _  i  )  I 


j  1 

I  (*k- An n  .e.'l. ,0 ~tc .  1  ■  •‘U  .0  ^ I 


where 


(8) 


26u7  (4,o) 

6  >  0 

57  (tfc,o)  €  R 


After  b  and  zu  (4,o)  are  found  from  the  above  system,  a  (tk, o)  is  computed  using 
equation  4. 

In  a  (tfc.o)  =  In  |  (sk,  ha k.0,tk,n,wk.,n,ak,0^  \  +  b(zu  (tk,o)  -  uk,m)2  (9) 

Similarly,  by  applying  the  test  points 

(«fc,0’  4,i>  k’fe.-i)  Pkx>) 

<*kfi,tkti,b*kfl,Pkfi)  (10) 

ftfc.O:  'h,o) 

another  set  of  intermediate  variables  w  (4,i)  and  a(tki, ,)  are  obtained.  Letting 

«fc,0/8fc  +  ,ln 

r  = - ; -  (11) 

<*fc  +  C*fc,0 

and  substituting  into  equation  4  gives 


W  (4,™)  =Uk-r  (4  -  4,») 


(12) 
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Using  w  (4, o)  and  w  (4,  1),  the  variable  r  can  be  computed. 


^  (4,o)  ^  (4,i) 

r  = - 

4,o  —  4,i 

Now,  the  parameters  a k  and  (3k  can  be  determined  from  the  following 


(13) 


b  = 


r  = 


_ ft  A: +<**-. 0 _ 

2  [(ftt +«*■.©)  +(j3k-i3ks>)  j 

nk&k-,0+Cck\0&k 

ftfc+ftfc,  0 


(14) 


Using  mathematical  manipulations,  the  following  is  obtained. 


2  b 


(0  k  0  k,  0) 


Oik,  0 

r  ~  £fc,0 


(/^fc  0k,  0)  —  0 


(15) 


The  solution  to  the  above  equation  is  given  by 


3  1 

ft|,0 

;.  =  < 

26(r-^.,0) 

'C-.JV 

(16) 


If  3k  =  3k  0,  then  from  equation  1 0.62  one  gets 


r  =  ■%,  0 


(17) 


and 


Otherwise, 


a“=2 b~  ° 


(18) 


3i.  —  r 

ak  ~  — ftfc.O 

T  Jk.  0 
>  0 


(19) 
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Now  that  ak  and  3k  are  known,  the  time  parameter  tk  can  be  computed  from  10.53 
as  follows. 


^  _  01  k 1  +  <4,o  ^  Q  (4,o)  _  4,0  +  4,1 
' k  4,o -4,i  na(4,i)  2 

Finally,  the  frequency  parameter  uik  can  be  solved  from  10.63 


(20) 


uik  =  w  (4,0)  +  r  (4  -  4.o)  (21) 

Recent  studies  show  that  radar  signal  can  be  properly  modeled  using  the  chirplet 
basis  function.  Many  natural  phenomena,  for  instance,  signals  encountered  in  radar 
systems  [26-28],  the  impulsive  signal  that  is  dispersed  by  the  ionosphere  [29],  and 
seismic  signals  [30],  are  modeled  as  chirplet  type  functions.  The  chirplet  function 
has  many  striking  characteristics.  1)  The  chirplet  basis  is  a  well  understood  four- 
parameter  function,  ak,0k,tk,  and  04,  localized  in  the  joint  time-frequency  plane. 
Only  a  moderate  computation  time  is  needed  to  search  for  the  basis  parameters. 
1)  Since  both  the  amplitude  modulation  (AM)  and  frequency  modulation  (FM)  are 
part  of  the  basis  function,  the  chirplet  can  more  efficiently  represent  the  reflected 
signal  from  a  target  with  a  rotating  part.  Only  a  few  set  of  these  bases  are  needed  to 
approximate  the  time-frequency  structure  of  the  radar  signal.  3)  More  importantly, 
the  radar  returns  from  the  target  body  and  the  rotating  part  can  be  more  easily 
separated  based  on  the  parameters  of  the  chirplet  bases,  particularly  nk  and  Ak. 
This  is  because  signals  from  the  main  body  and  the  rotating  part  are  captured  by 
chirplet  bases  with  different  parameters.  We  separate  the  body  signal  from  the 
rotating  part  using  threshold  values  of  a*,  and  3k.  After  the  separation,  we  process 
the  main  body  signal  and  the  rotating  part  signal  individually  for  better  information 
extraction.  This  includes  both  the  extraction  of  the  geometrical  features  from  the 
main  body  and  the  m-D  features  from  the  moving  parts. 

Following  the  extraction  of  m-D  features  from  the  radar  returned  signal  using  either 
wavelet  or  chirplet  analysis,  the  time-frequency  signatures  of  the  m-D  features  can 
be  utilized  to  visualize  the  oscillation  and  to  extract  motion  parameters  related  to 
the  target  of  interest.  The  reason  for  this  is  that  the  joint  time-frequency  analysis 
can  provide  time-varying  Doppler  frequency  information.  From  the  joint  time- 
frequency  representation  of  the  m-D  features,  the  motion  parameters  of  a  target 
such  as  the  period  and  amplitude  of  the  oscillation  can  be  estimated  [1 1-13]. 

To  further  the  analysis  of  the  extracted  m-D  features,  the  autocorrelation  of  the 
m-D  time  sequence  signals  can  be  computed  and  investigated.  The  autocorrela¬ 
tion  of  time  sequence  signals  has  already  been  used  by  Bell  and  Grubbs  to  es- 
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timate  the  vibration/rotation  rate  of  micro-motions  from  a  signal  containing  both 
vibrational/rotational  and  stationary  parts  [31].  It  is  probable  that  the  estimated 
rate  will  be  more  accurate  if  this  approach  is  used  on  data  including  only  vibra¬ 
tional/rotational  components  obtained  from  the  wavelet  decomposition.  This  is  be¬ 
cause  the  return  from  the  stationary  parts  strongly  overlaps  with  the  return  from 
the  vibrational/rotational  parts  in  the  Doppler  frequency  spectrum.  Namely,  signals 
from  stationary  parts  and  rotating  parts  superimpose  on  each  other  and  thus  are 
separable.  The  extracted  m-D  features  are  employed  to  estimate  the  target’s  motion 
parameters. 
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3  Micro-Doppler  Analysis 


The  sample  data  sets  are  APY-6  collected  at  PAX  for  a  rotating  antenna.  The  para¬ 
meters  of  the  SAR  system  used  to  collect  this  set  of  data  are  shown  in  Table  1.  The 
data  has  been  pulse  compressed  and  motion  compensated  to  the  first  level.  This  is 
an  initial  data  distribution  and  there  may  be  undiscovered  flaws  in  the  processing. 
A  Hamming  window  is  applied  to  suppress  range  side-lobes.  One  thousand  twenty 
four  samples  were  processed.  However,  the  data  cuts  may  not  contain  edge  samples 
that  are  reduced  by  the  IF  filter.  Full  data  cuts  are  72  seconds  long  and  there  are 
2048  pulses,  1024  range  bins  and  8  bytes  per  complex  sample.  Using  this  data  set 
and  employing  the  proposed  method,  the  rotation  rate  of  the  antenna  is  estimated. 
The  rotation  rate  of  the  rotating  antenna  is  known  to  be  4.7  seconds. 

The  original  SAR  image  is  shown  in  Figure  3.  The  Doppler  smearing  due  to  the  ro¬ 
tating  parts  is  often  well  localized  in  a  finite  number  of  range  cells.  It  is  reasonable 
to  process  the  Doppler  signal  for  each  range  cell  independently.  Since  the  ground 
truth  of  the  target  is  already  known,  the  data  between  the  range  cells  121  and  125 
was  analyzed  using  the  wavelet  decomposition  method  and  the  adaptive  chirplet 
decomposition  method.  Therefore,  the  discussion  of  this  report  will  deal  exclus¬ 
ively  with  these  range  cells.  Figure  4  shows  the  zoomed  in  SAR  image  between  the 
range  cells  1 15  and  130. 

The  following  two  subsections  give  the  results  obtained  after  signal  analysis  and 
processing  using  both  the  wavelet  decomposition  method  and  the  adaptive  chirplet 
decomposition  method  for  range  cells  123  and  124. 


3.1  Wavelet  Analysis 

The  analysis  of  the  Doppler  signals  using  the  wavelet  decomposition  method  is 
given  in  Figures  5  through  18  for  range  cells  123  and  124.  The  analysis  of  each 
range  cell  includes: 

1 .  Time-series  of  the  original  signal, 

2.  Fourier  transform  of  the  original  signal, 

Table  1:  The  characteristics  of  SAR 


Center  frequency 

9.65  GHz 

Pulse  repetition  Interval  (PRJ) 

1/1200  seconds 

Bandwidth 

150  MHz 
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Figure  3:  The  original  SAR  image. 

3.  Time-frequency  signature  of  the  original  signal, 

4.  Time-frequency  signature  of  the  extracted  body  signal, 

5.  Time-frequency  signature  of  the  extracted  oscillating  signal, 

6.  Autocorrelation  function  versus  time  for  the  extracted  body  signal,  and 

7.  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal. 

The  wavelet  transform  method  consists  of  decomposing  the  original  signal  into 
components  and  recombining  these  components  in  various  groupings  in  order  to 
obtain  the  reconstructed  signals  that  represent  the  desired  features.  In  this  case,  the 
features  extracted  are  the  body  signal  and  the  oscillating  signal. 
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Figure  4:  The  zoomed  in  SAR  image. 
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The  Time-Senes  of  the  Ortgmal  Data 
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Figure  5:  Time-series  of  the  original  signal  at  range  cell  123 
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Figure  6:  Fourier  transform  of  the  original  signal  at  range  cell  123 
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Spectrum  TF  Signature  of  the  Original  Signal 


Figure  7:  Time-frequency  signature  of  the  original  signal  at  range  cell  123 


16 


DRDC  Ottawa  TM  2005-204 


Figure  8:  Time-frequency  signature  of  the  extracted  body  signal  at  range  cell  1 23 
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Figure  9:  Time-frequency  signature  of  the  extracted  oscillating  signal  at  range  cell 
123 
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Figure  10:  Autocorrelation  function  versus  time  for  the  extracted  body  signal  at 
range  cell  123 


18 


DRDC  Ottawa  TM  2005-204 


The  Autocorrelation  of  the  Oscillating  Signal 
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Figure  11:  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal 
at  range  cell  123 


The  Time-Series  of  the  Original  Data 
1  - - 


09 

0.8 

0.7 


Time  (sec) 

Figure  12:  Time-series  of  the  original  signal  at  range  cell  124 
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The  Fourier  Spectrum  of  the  Original  Data 
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Figure  13:  Fourier  transform  of  the  original  signal  at  range  cell  124 


Spectrum 


TF  Signature  of  the  Original  Signal 


Figure  14:  Time-frequency  signature  of  the  original  signal  at  range  cell  124 
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Spectrum 


TF  Signature  of  the  Body 


Figure  15:  Time-frequency  signature  of  the  extracted  body  signal  at  range  cell  124 
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Figure  16:  Time-frequency  signature  of  the  extracted  oscillating  signal  at  range  cell 
124 
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The  Autocorrelation  of  the  Body  Signal 
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Figure  17:  Autocorrelation  function  versus  time  for  the  extracted  body  signal  at 
range  cell  124 
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Figure  18:  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal 
at  range  cell  124 
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3.2  Adaptive  Chirplet  Analysis 

The  analysis  of  the  Doppler  signals  using  the  adaptive  chirplet  decomposition  is 
given  in  Figures  19  through  32  for  range  cells  123  and  124.  Similar  to  wavelet 
analysis,  the  analysis  of  each  range  cell  includes: 

1.  Time-series  of  the  original  signal, 

2.  Fourier  transform  of  the  original  signal, 

3.  Time-frequency  signature  of  the  original  signal, 

4.  Time-frequency  signature  of  the  extracted  body  signal, 

5.  Time-frequency  signature  of  the  extracted  oscillating  signal, 

6.  Autocorrelation  function  versus  time  for  the  extracted  body  signal,  and 

7.  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal. 

As  with  the  wavelet  transform  method,  the  features  extracted  here  are  the  body 
signal  and  the  oscillating  signal. 
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The  Time-Senes  of  the  Original  Data 
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Figure  1 9:  Time-series  of  the  original  signal  at  range  cell  1 23 


The  Fourier  Spectrum  of  the  Original  Data 
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Figure  20:  Fourier  transform  of  the  original  signal  at  range  cell  123 
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Spectrum 


TF  Signature  of  the  Original  Signal 


Figure  21:  Time-frequency  signature  of  the  original  signal  at  range  cell  123 


Spectrum 


TF  Signature  of  the  Body 


Figure  22:  Time-frequency  signature  of  the  extracted  body  signal  at  range  cell  123 
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Figure  23:  Time-frequency  signature  of  the  extracted  oscillating  signal  at  range  cell 
123 


The  Autocorrelation  of  tha  Body  Signal 


Figure  24:  Autocorrelation  function  versus  time  for  the  extracted  body  signal  at 
range  cell  123 
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The  Autocorelation  of  the  Oscillating  Signal 


Figure  25:  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal 
at  range  cell  1 23 


Figure  26:  Time-series  of  the  original  signal  at  range  cell  124 
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The  Fourier  Spectrum  of  the  Original  Data 
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Figure  27:  Fourier  transform  of  the  original  signal  at  range  cell  124 


Spectrum  TF  Signature  of  the  Original  Signal 


Figure  28:  Time-frequency  signature  of  the  original  signal  at  range  cell  124 
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Spectrum 


TF  Signature  of  the  Body 


Figure  29:  Time-frequency  signature  of  the  extracted  body  signal  at  range  cell  124 


Spectrum  TF  Signature  of  the  Oscillation 


Figure  30:  Time-frequency  signature  of  the  extracted  oscillating  signal  at  range  cell 
124 
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Figure  31:  Autocorrelation  function  versus  time  for  the  extracted  body  signal  at 
range  cell  124 


The  Autocorrelation  of  the  Oscillating  Signai 


Figure  32:  Autocorrelation  function  versus  time  for  the  extracted  oscillating  signal 
at  range  cell  124 


30 


DRDC  Ottawa  TM  2005-204 


4  Analysis  Results 


The  results  obtained  from  both  the  analysis  using  the  wavelet  decomposition  method 
and  the  adaptive  chirplet  decomposition  method  allow  the  rotation  rate  of  the  an¬ 
tenna  to  be  estimated.  This  is  done  by  examining  the  time-series  of  the  original 
data,  the  time-frequency  signature  of  the  extracted  oscillating  signal,  and  the  auto¬ 
correlation  of  the  extracted  oscillating  signal  for  each  of  the  range  cells  (using  both 
methods). 

Using  the  time-frequency  plot,  the  rotation  rate  of  the  antenna  is  estimated  by  meas¬ 
uring  the  time  interval  between  peaks.  The  period  is  the  time  interval  between 
peaks.  As  an  example,  in  Figure  9,  there  are  3  peaks.  The  time  interval  between 
peak  1  and  peak  2,  between  peak  2  and  peak  3,  and  between  peak  1  and  peak  3 
are  measured.  The  average  value  is  then  used  to  estimate  the  rotation  rate  of  the 
antenna.  The  period  of  the  rotating  antenna  can  be  estimated  by  taking  the  auto¬ 
correlation  of  the  time  sequence.  The  peaks  in  correlation  coefficient  correspond  to 
the  dominant  period  of  the  signal. 

The  results  of  range  cells  121  and  125  are  not  clear  enough  to  provide  any  signific¬ 
ant  indications  towards  the  rotational  rate  for  either  method.  However,  for  the  data 
analyzed  herein,  the  rotational  rate  can  be  estimated  for  range  cells  122,  123,  and 
124.  These  values  are  given  in  Table  2. 

The  actual  antenna  rotational  rate  is  known  to  be  4.7  seconds.  Clearly,  the  estim¬ 
ated  values  for  both  the  wavelet  decomposition  method  and  the  adaptive  chirplet 
decomposition  method  for  range  cells  122,  123,  and  124  correspond  directly  with 
the  known  value. 

Overall,  range  cell  123  gives  the  best  results  for  both  methods.  This  is  evident  by 
looking  at  the  time-frequency  signatures  of  the  original  and  extracted  signals.  Spe¬ 
cifically,  the  oscillations  of  the  rotating  antenna  can  be  clearly  seen  by  looking  at 
the  time-frequency  signatures  of  the  extracted  oscillating  signal  in  Figures  9  and 
23  (the  wavelet  and  chirplet  methods  respectively).  Furthermore,  the  autocorrela¬ 
tions  of  the  extracted  oscillating  signal  in  Figures  1 1  and  25  (wavelet  and  chirplet 

Table  2:  Estimated  Rotation  Rate  of  the  Antenna  (seconds) 


Range  Cell 

Wavelet  Decomposition 

Adaptive  Chirplet  Decomposition 

122 

4.6 

4.4 

123 

4.8 

4.8 

124 

4.7 

4.8 
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respectively)  both  give  the  estimated  rotation  rate  to  be  4.8  seconds,  which  is  very 
close  to  the  actual  value  of  4.7  seconds. 
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5  Conclusion 


Moving  targets  cause  phase  modulation  of  the  SAR  phase  history  corresponding  to 
the  target.  A  special  case  is  rotating  point  targets.  This  report  examines  the  case 
of  a  rotating  antenna  within  in  a  SAR  target  scene  collected  by  the  US  Navy  APY- 
6  radar.  Rotating  target  causes  sinusoidal  phase  modulation  of  the  azimuth  phase 
history  of  the  SAR  collection.  The  phase  modulation  may  equivalently  be  seen  as 
a  time-varying  Doppler  frequency.  Time-frequency  methods  are  useful  tools  in  the 
analysis  of  such  signals  as  they  use  time  integration  while  still  allowing  for  non¬ 
stationary  signals. 

In  this  report,  we  present  two  approaches  for  m-D  analysis  for  the  extraction  of 
the  m-D  features  of  radar  returned  signals  from  targets  in  a  SAR  scene.  Both  ap¬ 
proaches  have  two  main  processes;  the  first  being  the  decomposition  of  the  SAR 
return  in  order  to  extract  m-D  features,  and  the  second  being  the  time-  frequency 
analysis  in  order  to  estimate  motion  parameters  of  the  target  in  the  SAR  image.  Two 
different  approaches  are  used  in  the  decomposition  process.  The  first  approach,  a 
wavelet  decomposition  analysis,  has  been  used  in  our  previous  m-D  analysis  work 
for  a  ISAR  system,  and  it  is  presented  here  for  the  extraction  of  m-D  features  for  a 
SAR  system.  The  second  approach,  an  alternative  to  the  wavelet  approach,  is  a  new 
approach  that  is  introduced  for  decomposing  the  time  domain  radar  signal  in  SAR 
based  on  the  adaptive  chirplet  decomposition.  Different  time-frequency  techniques 
are  employed  depending  on  the  procedure  used  in  the  decomposition  process.  For 
the  wavelet  based  approach,  the  short-time  Fourier  transform  is  utilized,  while  in 
the  procedure  based  on  the  adaptive  chirplet  decomposition,  the  time-frequency 
distribution  series  is  used. 

By  applying  the  proposed  procedures  to  data  from  a  rotating  antenna,  the  effect¬ 
iveness  of  both  these  analysis  techniques  in  a  SAR  target  scenarios  is  confirmed. 
From  the  extracted  m-D  signatures,  information  about  the  target’s  micro-motion 
dynamics,  such  as  rotation  rate,  has  been  obtained. 

As  demonstrated,  both  wavelet-based  and  chirplet-based  procedures  are  successful 
for  the  extraction  of  m-D  features  from  SAR  data.  The  wavelet  based  method  re¬ 
quires  the  decomposed  signal’s  constituent  parts  to  be  combined  in  a  specific  man¬ 
ner  in  order  to  extract  m-D  features.  While  this  is  not  the  case  for  the  chirplet  based 
method,  the  chirplet  based  method  still  requires  its  alpha  and  beta  parameters  to 
be  specifically  set  in  order  to  extract  any  m-D  features.  The  chirplet  based  method 
establishes  itself  to  be  a  viable  alternative.  In  general,  however,  it  is  shown  that  the 
results  are  much  improved  after  the  m-D  extraction  from  the  SAR  target  scenarios 
has  taken  place  since  only  the  vibration/rotational  components  are  employed.  The 
rotating  frequency  of  the  antenna  was  calculated  using  these  two  approaches  and 
the  results  agree  well  with  available  ground  truth.  Such  motion  parameters  may  be 
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useful  for  describing  the  target,  or  serve  as  inputs  to  automatic  target  recognition 
algorithms. 
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indexing  terms  which  are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 
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